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We present an analytical study of the D'yakonov-Perel' spin relaxation time for degenerate elec- 
trons in a photo-excited electron-hole liquid in intrinsic semiconductors exhibiting a spin-split 
band structure. The D'yakonov-Perel' spin relaxation of electrons in these materials is controlled 
by electron-hole scattering, with small corrections from electron-electron scattering and virtually 
none from electron-impurity scattering. We derive simple expressions (one-dimensional and two- 
dimensional integrals respectively) for the effective electron-hole and electron-electron scattering 
rates which enter the spin relaxation time calculation. The electron-hole scattering rate is found to 
be comparable to the scattering rates from impurities in the electron liquid - a common model for 
n-type doped semiconductors. As the density of electron-hole pairs decreases (within the degener- 
ate regime), a strong enhancement of the scattering rates and a corresponding slowing down of spin 
relaxation is predicted due to exchange and correlation effects in the electron-hole liquid. In the 
opposite limit of high density, the original D'yakonov-Perel' model fails due to decreasing scattering 
rates and is eventually superseded by free precession of individual quasiparticle spins. 



I. INTRODUCTION 

Semiconductor spintronic devices are poised to com- 
plement existing electronic devices in the future by using 
the electron's spin degree of freedom to store and transfer 
information. 1 ^ 2 Of the many branches and uses of spin- 
tronics, spin relaxation is relevant to each. The rate of 
decay of spin polarization needs to be taken into con- 
sideration in spin packet transportation, spin transistor 
junctions, spin qubits, and the like. Additionally, micro- 
scopic studies of spin relaxation severely test our ability 
quantitatively describe effective interactions and scatter- 
ing mechanisms in the solid state. 3-7 

The most complete study of spin relaxation in semi- 
conductors to date is by Jiang et al., 8 who made calcula- 
tions of the spin relaxation time (SRT) due to each of the 
relevant spin relaxation mechanisms in III-V semicon- 
ductors: Bir-Aronov-Pikus (BAP) 9 , D'yakonov-Perel' 
(DP) 10 , and Elliott- Yafet (EY) 11 . Comparing the rel- 
ative efficiencies of the mechanisms, their results suggest 
that D'yakonov-Perel' spin relaxation is dominant for es- 
sentially all the electron densities and temperatures of 
experimental interest. 8 Although Teng et al. initially 
found cases where BAP spin relaxation dominates DP 
spin relaxation in intrinsic GaAs, 12 it was later pointed 
out by Jiang and Wu that non-degenerate statistics were 
being applied to degenerate electrons. 13 Also noted by 
Jiang et al., Song and Kim have investigated spin relax- 
ation due to all of the relevant spin relaxation mecha- 
nisms in n- and p-type III-V semiconductors, 12 but their 
analytical expressions are again only applicable to the 
non-degenerate regime. Tamborenea et al. 14 have calcu- 
lated EY SRTs due to electron-impurity and electron- 
electron collisions in n-doped GaAs for a wide range of 
temperatures and electron densities, finding that the EY 



mechanism alone is insufficient to explain experimental 
SRT measurements, 15 suggesting that DP spin relaxation 
may account for the discrepancies. These findings en- 
courage us to continue looking at DP spin relaxation as 
the primary spin relaxation mechanism in GaAs and sim- 
ilar semiconductors of the zincblende structure. 

The focus of the present study is on the role played 
by many-body interactions on the DP spin relaxation 
mechanism - a role that we would like to clearly dis- 
entangle from that of other scattering mechanisms. The 
theory of many-body effects in DP spin relaxation for 
electrons in semiconductors was formulated by Glazov 
and Ivchenko (GI) several years ago. 16 ' 17 The basic idea 
is that electron-electron scattering causes the electron 
quasiparticle to perform a random walk in momentum 
space; this in turn causes random variations of the di- 
rection and magnitude of the spin precession axis. DP 
spin relaxation arises from the cumulative effect of many 
small precessions about randomly varying axes, and its 
main signature is that the spin relaxation rate is inversely 
proportional to the momentum scattering rate. In prac- 
tice, however, it is very difficult to disentangle the con- 
tribution of the electron-electron scattering rate from the 
similar but much larger contribution of electron-impurity 
scattering, which is inevitably present in doped semicon- 
ductors. An important exception arises in intrinsic semi- 
conductors, when a non-equilibrium population of elec- 
trons and holes can be created by optical excitation. By 
using circularly polarized light it is possible to achieve a 
high degree of spin polarization of the electrons in the 
conduction band (the spins of the holes relax rapidly due 
to strong spin-orbit interactions in the valence band), 
and the time evolution of this spin polarization can be 
monitored in real time. Such a system is virtually impu- 
rity free, and thus offers a unique opportunity to directly 
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test the impact of many-body interactions on DP spin 
relaxation of electrons. 

The initial motivation for the present study came from 
an experiment in which a density-dependent diffusion 
rate for spin polarized electrons was observed in a photo- 
excited electron- hole packet in bulk GaAs. 18 Zhao et al. 
found that the spin density in such a packet has a smaller 
diffusion constant than the carrier density. It was ar- 
gued that this could be explained by a density dependent 
spin relaxation rate, where the electron spins in low den- 
sity regions of the packet relax faster than in the high 
density regions, leading to the appearance of slow diffu- 
sion. This behavior is consistent with that of DP spin re- 
laxation controlled by electron-electron scattering in the 
non-degenerate regime, 17 which is indeed expected to oc- 
cur in the low density regions of the electron-hole packet. 
Remarkably, the increase of the spin-relaxation rate with 
decreasing density in the non-degenerate regime is oppo- 
site to the behavior in the degenerate regime, 8 where the 
spin relaxation rate decreases with decreasing density. 
This means that the SRT has a maximum as a function 
of density at a density intermediate between the degen- 
erate and the non-degenerate regime. 

The present paper continues our efforts to develop a 
better microscopic understanding of the GI mechanism 
of DP spin relaxation. We consider electron-electron, 
electron-impurity, and electron-hole interactions in the 
degenerate regime. In this regime, the presence of a 
large Fermi surface paves the way to an elegant analytical 
treatment of electron-electron and electron-hole scatter- 
ing rates - a treatment pioneered by Abrikosov and Kha- 
latnikov in their classic paper on the theory of the Fermi 
liquid. The relevant scattering rates are similar to but 
not identical with the well-known momentum scattering 
rates which control, for example, the lifetime of quasi- 
particles. The difference arises because different collision 
events are not equally effective at changing the magni- 
tude and direction of the spin precession: the contribu- 
tion of different collision events must be weighted accord- 
ing to their effectiveness at changing the spin precession. 
(A similar "weighting" - this time concerning the effec- 
tiveness of collisions at changing the current - is responsi- 
ble for the replacement of the scattering cross section by 
the "transport cross section" in microscopic calculations 
of the electrical conductivity). GI had previously looked 
at the contribution of electron-electron collisions to DP 
spin relaxation in the degenerate regime, but limited 
themselves to qualitative estimates. 17 Jiang, et al. ap- 
plied the powerful kinetic spin Bloch equation approach, 
including all relevant spin relaxation mechanisms, but 
their work was primarily numerical. 8 In the present pa- 
per, analytic expressions are derived which are appli- 
cable to a variety of semiconductors of the zincblende 
structure; in particular, the calculation of the relevant 
electron-electron and electron-hole scattering rates is re- 
duced to the evaluation of simple two-dimensional and 
one-dimensional integrals over the Fermi surface. We find 
that electron-hole scattering is the dominant mechanism 



of momentum randomization in these intrinsic semicon- 
ductors. Nevertheless, we emphasize that our mechanism 
of spin relaxation remains conceptually distinct from the 
BAP mechanism in which electrons transfer their spin 
polarization to the holes via interband matrix elements 
of the Coulomb interaction. In the present mechanism 
- that is in the DP mechanism - electrons transfer only 
momentum to the holes. 

Another novel feature of the present work is that we 
include exchange and correlation effects in the calcula- 
tion of the effective electron-electron and electron-hole 
interactions. In practice, this means that we are going 
beyond the basic RPA-screened interaction. Exchange 
and correlation effects arc included via local field factors, 
and we will show that these effects cause an enhancement 
of the scattering rate - and a corresponding reduction of 
spin relaxation rate - at low density and low temperature 
(still in the degenerate regime). In the opposite limit of 
high density, the scattering lifetime becomes very large. 
A cautionary note is included for the treatment of DP 
spin relaxation in this regime. Specifically, the common 
assumption that momentum-changing collisions are fre- 
quent on the scale of spin relaxation fails, since the two 
processes occur on comparable time scales. In this regime 
the quasiparticles are essentially non-interacting and the 
DP mechanism is superseded by the spin precession of 
individual quasiparticles of essentially constant momen- 
tum. 

This paper is organized as follows: in Section II, we 
review the steps and the assumptions in the derivation 
of the standard formula for D'yakonov-Perel' spin 
relaxation; in Section III we derive analytic expressions 
for the effective (i.e. weighted) scattering rates due to 
electron collisions with electrons, holes and impurities: 
Section IV describes the effective scattering amplitudes 
used in our calculations; Section V discusses the effective 
scattering rates and SRTs for doped and intrinsic GaAs; 
Section VI contains our concluding remarks. 



II. D'YAKONOV-PEREL' SPIN RELAXATION 

Although the derivation of DP spin relaxation is well 
known and can be found in excellent reviews 3,5 we briefly 
reproduce it here, partly to adapt the notations, and 
partly to spell out the underlying physical assumptions. 
The system under investigation is a III-V semiconduc- 
tor (i.e. GaAs, InAs, InSb) of the zibcblende structure, 
which exhibits spin-split bands. The spin splitting is 
caused by a Drcssclhaus effective magnetic field 19 which 
arises from spin-orbit interaction and couples to electrons 
via the hamiltonian term 

^ik-^k-tr, (1) 

where £1^ defines the precession axis and precession fre- 
quency of electron spins. f2k is a cubic harmonic function 
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of the Bloch wave vector k; its component along the z- 
axis is given by 



k z (k x k y ) , 



(2) 



where a c is the Dresselhaus spin-orbit coupling constant 
(e.g. a c = 0.07 for GaAs), E g is the band gap energy, 
and m c is the conduction band effective mass of an elec- 
tron. The other components of i7 k are obtained by cyclic 
permutations of x, y, and z. 

The spin Boltzmann kinetic equation describes the evo- 
lution of electrons in time: 



dp k 1 dp k 



F k +^-v k = / k (i), (3) 



where p k is the 2x2 spin density matrix and i k (t) is 
the collision integral (also a 2 x 2 matrix). Assuming a 
homogeneous distribution of electrons (dp^/dr — 0) and 
absence of external fields (F k = 0), Eq. (3) reduces to 



dpk 1 [n 



o",Pk] = Iu{t) ■ 



(4) 



We prepare a spin polarized distribution which evolves 
according to Eq. (4) and eventually relaxes to a com- 
pletely unpolarized state. To describe the process, we 
search for a solution which contains a quasi-equilibrium 
component, p ok , describing a state of uniform spin po- 
larization along a direction denoted by s, and a non- 
equilibrium component, pi k : 



Pk = Pok + Pik ■ 



(5) 



Using this definition in Eq. (4), we can relate terms of 
varying order of spin orbit interaction (a c ): 



dp, 



'Ok 



dt 



+ 



dt 2i l 



! POk] 

- ^[flk • er,Pik] = Ik(t) ■ 



(6) 



The following presumptions are made about the various 
terms appearing in Eq. (6): 



1. Order of spin-orbit interaction (a c ): 
th order; p lk : 1 st order; 
2 nd order; pi k : 3 rd order. 



Pok- 

POk- 



2. I k (i) may be cast in the form of relaxation time 
approximation. 

The first of these presumptions will be shown to be con- 
sistent with the form of the kinetic equation momentar- 
ily. The legitimacy of the second presumption will be 
discussed in the next section when we analyze the colli- 
sion integral in detail. For now, we simply write 



Pik 



(7) 



where t£ is an effective scattering time on the order of 
the plane wave lifetime. The precise form of depends 
on the scattering mechanism and will be explicitly con- 
structed in the next sections for various cases of interest. 
Notice that pok, being a quasi-equilibrium distribution, 
does not contribute to the collision integral. 

Equating terms of the same order in a c in Eq. (6), the 
following relations emerge: 



dpok 

dt 
dpik 

dt 



Pik 



[fi k • <T,P0k] 



2i 
0. 



7 ([fJ k • er,p lk ]) fc 



(8) 
(9) 
(10) 



where (■ ■ - )k stands for the angular average over all di- 
rections of k at fixed magnitude k — |k|. The angular 
averaging is necessary in Eq. (9) to properly relate the 
quantities. Finally, dpik/dt has no counterpart in the 
kinetic equation at order a^, hence it must vanish at this 
order as anticipated. 

Let us now introduce the form of the quasi-equilibrium 
distribution, po k , from which we may use the relations 
above to find expressions for p lk and the spin polariza- 
tion, pofe is written in terms of equilibrium Fermi-Dirac 
distributions for spin-up electrons (f') and spin-down 
electrons (/■'■) in the presence of a fictitious fc-dependent 
magnetic field which we introduce to enforce the desired 
quasi-equilibrium spin polarization: 



Pok = fk 1 + /feS • c 
where 1 is the 2x2 identity matrix, 

m >' k 1 >> k 

.Ik o 



(11) 



(12) 



is the spin-averaged occupation number of the k state 
and 



£S _ J k J k 

Jk 2 



(13) 



is the quasi-equilibrium spin polarization occupation 
number of the k state. The distributions for spin-up and 
spin-down electrons are defined by 



Jk 



i 



1 -f e l3{e k -^e s ) 



(14) 



where j3 = l/(kgT), /i is the chemical potential, and e s 
is the Zeeman interaction energy. 

The non-equilibrium part of the distribution function 
arises from the competition between spin precession, 
which tilts the spins away from the s-axis, and collisions, 
which attempt to restore local equilibrium. Inserting pofe 
[Eq. (11)] into the relation for pi k [Eq. (8)], we find 



Pik = (ft k x /|s) • a . 



(15) 
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At this point, the presumption made earlier in this sec- 
tion concerning the orders of a c in pok, Pik, Pofc, and pik 
may be verified easily. 

We identify the quasi-equilibrium spin polarization at 
wave vector k as 



and notice that it does not depend on the direction of k. 
Then, the time rate of change of s fc is found by tracing 
Eq. (9) with er: 



9s, 



Ofc 



dt 



(Tk^k X (O k X S fc))fc . 



(17) 



This equation can be cast into the standard relaxation- 
type form, 



9s, 



Ofc 



dt 



s 0k 



(18) 



by defining the rate of spin relaxation at k as follows: 
1 



» 



(r^-^k-s) 2 )),. 



(19) 



The cubic symmetry of III-V semiconductors implies 
that (i) the relaxation rate of Eq. (19) is independent of 
the direction of s, and (ii) the effective collision time 
depends on k only through its modulus k (in the next 
sections we verify this explicitly). Therefore, for cubic 
systems, the expression for the spin relaxation rate (19) 
simplifies further 



(20) 



The experimentally measurable quantity of primary 
physical interest is the expectation value of the spin den- 
sity 



s = s ofc . 

k 



(21) 



The equation of motion for S(t) is obtained by summing 
Eq. (17) over all k and using this result in the relaxation 
time approximation defining t^: 



OS 

~dt 



S 



(22) 



where the physical spin relaxation time is calculated from 
1 _2Ek«>* 



r ( 8 ) 



Ek fk 



(23) 



Obviously, the spin relaxation time as an internal 
material parameter, independent of a particular non- 
equilibrium state, is only meaningful in the regime of 



small spin polarization. In this limit, we can linearize /| 
in the polarization energy, s s : 



fk 



dft 
' de k 



(24) 



where ft is the unpolarized Fermi-Dirac function. Sub- 
stituting Eq. (24) into Eq. (23) we indeed find that the 



polarization energy, e s , cancels and we are left with the 
following expression which contains only internal charac- 
teristics of the system: 



rW 3 Ek/ fe °(l-/ fe °) 



(25) 



The spin relaxation time defined by this equation is 
the experimentally accessible quantity which determines 
a universal long time tail, S(t) ~ e _t / T< , in the spin 
relaxation dynamics. 

In the degenerate limit, k B T <C Sp, the factor ft(l — 
ft) entering the intergals in Eq. (25) is strongly peaked 
at the Fermi level. Since the spin orbit field ilk in Eq. (2) 
is a slowly varying function of k, it is legitimate to set 
k = kF and take it out of the integral. As a result, the 
spin relaxation rate takes the following simple form 



1 



g r avg 



where 



Ek T fc/fc(l ft 



avg 



Ek/ fe °(i-/ fe °) 



(26) 



(27) 



Hence, in the degenerate limit, the spin relaxation time 
is completely determined by two parameters: the 
mean square of the spin-orbit field (^}k F at the Fermi 
energy and a properly averaged effective scattering time 
r* vg . While calculation of (£l k ) fef , is straightforward, to 
find r* vg we normally need to solve a complicated integral 
equation. One of the goals of this paper is to show that 
r* vg , and thus , can be calculated rigorously using 
methods developed in the theory of Fermi liquids. 20 

There are two important points to be made about 
the above derivation. First, the whole treatment is re- 
liant on the spin-orbit interaction being weak, leading 
to a clear separation of time scales between the mi- 
croscopic momentum changing-collisions (fast) and the 
macroscopic spin relaxation (slow). In other words, we 
must have f2kTk <C 1; spins relax due to a combination 
of k-randomizing collisions and relatively small preces- 
sion rotations. From here, we can see that when the 
scattering lifetime is of the same order of magnitude as 
the spin-relaxation time, our treatment is no longer jus- 
tified. If spins relax just as quickly as collisions occur, 
it makes no sense to introduce a quasi-equilibrium dis- 
tribution function. In the extreme limit of infrequent 
collisions, the momentum of the electron becomes a con- 
stant of the motion and the spin simply precesses in the 
Dresselhaus field at a given k. 
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The second point is that Eq. (19) is often estimated by 
replacing the effective scattering time with the plane 
wave lifetime Tk- This, however, is not correct, for the 
reasons described in the introduction. r£ contains an ad- 
ditional momentum axis re-orientation mechanism, due 
to the Dresselhaus field. In the next section we will ob- 
tain the correct expression for t£, and in Section V we 
will explore the validity of the approximate replacement 
of t* by r k . 

III. THE EFFECTIVE SCATTERING TIME 

In order to calculate the effective scattering time rjj: 
we need to construct the collision integral. Naturally, 
the construction depends on the nature of the collision 
process. Let us begin with electron-electron scattering 
processes. Electron-hole scattering processes will then 
be easily handled. These two scattering mechanisms arc 
all that is needed to treat DP spin relaxation in intrinsic, 
photo-excited semiconductors. For doped semiconduc- 
tors, electron-impurity scattering needs to be included. 

A. Electron-electron collisions 

The collision integral for an interacting Fermi liquid 
has been derived by several authors. 17,21 ' 22 The most di- 
rect derivation starts from the Kadanoff-Baym quantum 
kinetic equations. 23 The collision integral derived in this 
manner has the form 

h(t) = 4 n (i) - I?* {t) , (28) 

where 

J — oo 

+G>(t,t 1 )X<(t 1 ,t)} , (29a) 
J£ ut (f)= / dti {££(*, ti)G£(ti,t) 

J — oo 

+G<(t,h)X>(t u t)} . (29b) 

Here, the lesser and greater Green's functions and the 
corresponding self-energies are all 2 x 2 matrices in spin 
space. For the self-energy function one adopts the Born 



approximation, in other words, second order in an effec- 
tive electron-electron interaction. The scattering ampli- 
tude is denoted by W'kpk'p' , where k, p are the momenta 
of the incoming particles and k',p' those of the outgo- 
ing ones (Notice that one must have k + p = k' + p' by 
momentum conservation). In the simplest approxima- 
tion W is simply the Fourier transform of the Coulomb 
interaction: 

47re 2 

Wkpk'p' = e | k _ k ,| 2 frc+p,k'+ P ' . (30) 

where e is the background dielectric constant of the semi- 
conductor and e is the electron charge. Here, we will in- 
clude RPA screening and post-RPA screening via local 
field factors. In a more complete treatment, the effective 
interaction would also depend on the relative orientation 
of the spins; we will not consider such effects here, but 
simply assume that W is averaged over the relative spin 
orientations. 

In Fig. 1 we show the diagrams that contribute to the 
collision integral in this approximation. As a further 
approximation we adopt the generalized Kadanoff-Baym 
(GKB) ansatz 23,24 relating the Green's functions to the 
distribution function: 

G£(t, = 0k(t)G£(t, - G£(t, i'K(i') , (31) 
where 

G r k {t, t 1 ) = -i0(t - t') e -^(t-t') ; ( 32a ) 

G£(t, t') = i6(t' - ty^-*") , (32b) 

and = h 2 k 2 /2m* is the single-electron energy in a 
parabolic band of effective mass m* . Notice that we are 
neglecting interaction contributions to the single particle 
energy, coming, for example, from exchange. It is as- 
sumed that these contributions are properly included in 
the effective mass. As a final approximation we ignore 
the difference between Pk(ii) and pk{t), i-e. we assume 
that the distribution function is slowly varying on the 
time scale probed by the integrals in Eq. (29). This is 
commonly referred to as the Markovian approximation, 
and is completely justified in the solution of a steady- 
state problem. Then, the remaining integral over time 
contains only a principle value term and a 5 of conser- 
vation of energy. Retaining only the dissipative contri- 
bution (from the conservation of energy term), the final 
expression is 



Ik{t) C ~ C — — — ^ 6 (Sk + £p — £k> — £p') ^k+p.k'+p' 
k'pp' 

[l^kpk'pf ({Pk, 1 - Pk'} Tr [p p (i - p p ,)] - {i - p k ,Pk'} Tr [(i - p p ) p p ,]) 

- iVFkpk'p'W^pp'k'l {{pk, (i - Pk') Pp (i - p P ')} - {i - Pk,/?k' (i - Pp) Pp'})] ■ (33) 



FIG. 1. (Color online) Contributions to the collision integral under the Born approximation. 



We now show that the relaxation-time approximation 
in Eq. (7) is consistent with this collision integral lin- 
earized in puc. With S (e^ + e p — £y — e p /), the quasi- 
equilibrium distributions make no contribution: 



{p ok , (l - pofc')} Tr [Po P (1 - Pop')] = 

{ (i - Pofc) , Pofc' } Tr [(1 - p 0p ) Pop'] (34) 



and 



and 



{Pok, (1 - Pok') Pop (1 - Pop')} = 

{(i - Pok) , Pofc' (i - Pop) Pop'} ■ (35) 



Then, each term in the collision integral contains one p\ : 



{pk, i - Pk'} Tr [p p (l - p p/ )] 
{1 - Pk,Pk'} Tr [(1 - p p ) p p /] = 

{pik, i - Pok'} Tr [p p (1 - Pop')] 
+ {pik,Pofc'}Tr [(1 - pop) pop'] 

- {pofc, Pik'} Tr [p p (i - Pop')] 

- {l - Pofc, Pik'} Tr [(1 - pop) Pop'] 
+ {pok, 1 - Pofc'} Tr [pi p (l - pop/)] 
+ {i ~ Pok,Pok'}Tr[pi P Pop'} 

- {pofc, 1 - Pofc' } Tr [popPip'] 

- {1 - p fe,Pofc'} Tr [(i - pop) pip'] 



(36) 



(37) 



{pk, (i - Pk') p p (i - p P ')} 
{i - pk,pk' (i - p P ) Pp'} = 

{pik, (i - P0k') Pop (1 - POp') } 
+ {pik, POfe' (i - POp) POp' } 

- {P0fe,Plk'P0p (l - POp')} 

- {i - Pofc,pik' (i - pop) pop'} 
+ {pok, (i - pot') pip (i - pop')} 

+ {l - P0fc,P0fc'PlpP0p'} 

- {POfe, (l - POfc') POpPlp'} 

- { i - Pofc, po/c' (i - pop) Pip' } • 



These can be simplified a bit further, but at a loss of 
readability. 

We assumed in Eq. (8) that the collision integral is 
proportional to pik, which is in turn proportional to 
[fit • er, pofc]- Furthermore, each component of J7k gets 
all its angular dependencies from / = 3 spherical harmon- 
ics [e.g. fik,z is expressed this way in Eq. (42)]. Obviously 
Pofc contains no such terms. Therefore pik must contain 
only / = 3 spherical harmonics. Viewing the collision 
integral as an integral operator, 



£pik = ^2^(K^)Pik' , 



(38) 



it is apparent that C, being rotationally invariant, has no 
means of changing the harmonic content of pik' • Then, 
we reach the conclusion that the collision integral also 
contains only / = 3 harmonics. In other words, we 
see that the relaxation time approximation introduced 
in Eq. (7) is consistent with the form of the linearized 
collision integral. 

Clearly, pik is traceless; so then is the collision integral. 
We trace the linearized collision integral with a to find a 
calculable relationship between elements in pik and Ik- 
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In the limit of small spin polarization p k is expanded to 
first order in the polarization energy with the assumption 
that e s -C e k . pok and pi k then take the forms 

Pak = f k l + Pe s f k \l- f k )s-cr, (39) 
Pik - K/3e s f k (l - f k ) (n k xs).<r, (40) 



where is the unpolarized Fermi-Dirac distribution. 
Equating terms of first order in e 8 we find 



/ fe °(l - /°)fi k x s = ^ ^ £ | dk' dp / fe °(l - /° - /° )<S (e k +e p - e k , - V ) <W' +P ' 

jl^kpk'pf bin k - r;,n k /] - \ |w kpk /p/W k p P / k< | [7-;n k - r;,n k / + T ;n p - T ;,n p ,] j x g. (4i) 



Cubic symmetry in III-V semiconductors permits us 
to consider only one component of this vector relation, 
f2 k x s — > O k z . The later quantity can be conveniently 
written in terms of spherical harmonics, Y l m (9, (j)): 



A entering Eq. (43) are defined as follows: 



„ft 2 fc 3 



k,z 



v / 2mp 9 



u)5 [Y£(# k ,<p k )+Yf 2 (0 k ,<p k j\ 



(42) 

This fact, together with isotropy of the scattering, imply 
that the solution of Eq. (41) does not depend on the 
direction of k. Indeed, assuming = r k , we find that 
the presence of fi kjZ [given by Eq. (42)] in the integrals 
in Eq. (41) already guarantees that the whole right hand 
side has the proper angular dependences, consistent with 
that of the left hand side. As further evidence of this 
fact, it is easily demonstrated that integration of f2 qjZ 
(where q can be any of k', p, or p') over c?q results in a 
term proportional to f2 k ,z- This is shown in Appendix A. 



As usual, further simplifications come in the degener- 
ate limit because the factor f k (l — f k /)fp(l — fpi) con- 
fines the momentum integrals to a narrow shell around 
the Fermi energy, where the density of states can be well 
approximated by a constant. In this case, T k = t*(£) 
becomes a function of the dimensionless energy variable, 
£ = (e k — /u)/(fcsT), which satisfies a one-dimensional 
integral equation of the following form: 



B 



f 

J — C 



(43) 



The first term of the integral corresponds to the com- 
bination of each of the first terms in square brackets of 
Eq. (41); the second term in the integral picks up all of 
the remaining terms of Eq. (41). A detailed derivation 
of this equation, starting with Eq. (41), is given in Ap- 
pendix B. The kernel K{x,^) and the parameters B and 



j £ ~ x 

7°(-0 U-e*-« 



= A 1 X 1 -A 2 X 2 /2 
A 1 - A 2 /2 ' 

J cos(0/2) 

J cos(fl/2) ' 

A 2 J cos( 



-A 2 J- cos(fl/2) [P3(C ° S01) 
-P 3 (cos9) + P 3 (cos9 2 )} . 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 

(50) 



P 3 (x) is the 3rd order Legendre polynomial. The solid 
angle <Kl = sm9d9d(f) should not be confused with 
the Dresselhaus Larmor frequency. All of these angles 
(6, (j), #i,#2) are described in Appendix B, but briefly: 
9 is the angle between k and p; <f> is the polar an- 
gle about the k + p axis, between k' (p') and k (p); 
cos#i = k • k' = (1/2)(1 + cos# + cos0 — cos#cos</>); 
cos 2 = k • p' = (1/2)(1 + cos 9 — coscj) + cos 9 cos <j>). 

Integral equations of the type in Eq. (43) with a ker- 
nel of Eq. (44) are common in the theory of transport 
coefficients of Fermi liquids . 25 » 26 The general method 
of solution has been proposed by Sykes and Brooker. 20 
This amounts to a clever Fourier transform utilizing the 
convolution theorem and ultimately changing the inte- 
gral equation to a recognizable inhomogeneous differen- 
tial equation. The solution to our Eq. (43) may then 
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literally be read out from their paper: 

cosh(£/2) f°° ^ _ iMi B^ (4Z + 3)<M") 

A 2( (A 2; -A) ' 
(51) 



r*(0 
where 



^ J-oo * & 



= pj +1 (tanhTrw) , 
A, = ^(i + l)(i + 2), 



(52) 
(53) 



and p™(x) are the associated Legendre polynomials. 

In Section II we showed that the physical spin relax- 
ation time is proportional to an average effective scatter- 
ing time r* vg given by Eq. (27). Expressing the summa- 
tion over k in Eq. (27) in terms of a ^-integral and using 
the solution of Eq. (51), we arrive at our final result: 



B 



E 



41 + 3 



(54) 



2tt 2 f - A 2l [A 2l - A] 



As it turns out, for all calculations performed in this 
paper, the first two terms of the sum account for greater 
than 99% of total scattering time. The first term alone 
is about 95% of the sum at the lowest density examined 
(n <~ 10 16 cm~ 3 ) and the accuracy increases with increas- 
ing density. 

For reference, we also mention the plane wave scatter- 
ing time. This can be calculated readily from Eq. (B19) 



by setting Pz{x) 
scattering time is 



and tZ — > r k - The fc-dcpcndcnt 



r fe = 25 



£fc - A< 



(55) 



and the averaged scattering time is 



Tavg ~ 2tt 2 



(56) 



Notice the close relationship between the plane wave scat- 
tering time and the I = term of the effective scattering 
time: 



B. Electron- hole collisions 

The collision integral for electron-hole collisions is sim- 
ply the direct-process-only portion of Eq. (33): 



71 

h 



-Z ^2 s ( £ k + £ P - £k> - v) <^k+p,k'+ P ' |Wkk'| 2 



k'pp' 

({p k ,i-p k '}Tr[pW (i-pW); 
- {i-p k ,Pk'}Tv[(i-pW)pW]) , (58) 

where is the hole density matrix. We consider equal 
densities of electrons and holes so that the Fermi mo- 
mentum wavevectors are equal for the two species. This 
condition is likened to an intrinsic semiconductor under 
optical excitation. Additionally, the valence band energy 
for holes is considered parabolic. 

The same logic used in solving the collision integral 
for electron-electron collisions (in Appendix B) can be 
applied here, with the modification 



dk' dp 



m c my 



2fi 6 cos(6»/2) 



dek' de p de p > sin 9d9 dip d<p p , 



where m v is the effective valence band hole mass. Fol- 
lowing the same procedure used in e-e scattering, the 
average effective scattering time for electron-hole colli- 
sions is 



4/ + 3 



2tt 2 ^ A 2l [A 2 , - Ai] 



(60) 



where the defines are the same as in Eqs. (45)-(50) and 

fi 7 (27T) 4 



B {h) = 



N 2 ^1 



m c ml (k B Tf 
The plane wave scattering time is simply 



with an average of 



£fc - M 
k B T 



r (h) 
avg 



3gW 
2tt 2 ' 



(61) 



(62) 



(63) 



'avg,Z=0 



w 

2^/r 2 



1 



1-A 



(57) 



For < A < 1, the effective scattering time is enhanced 
(or the rate is decreased) compared to the plane wave 
scattering time. This range of A is consistent with the 
calculations we make including the Dresselhaus field in 
Section V. 



C. Electron-impurity collisions 

The collision integral for elastic collisions from impu- 
rities is considerably simpler: 

/k(tr = | WW | 2 6(e k - e k ,) (p lk - pikO • 

n k' 

(64) 
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The resulting effective scattering rate is 3 ' 5 
1 ra c kFTii 

k F 



J d{co&e) \W{cos6)\ 2 [1 - P 3 (cos( 



where m is the density of impurities. 



(65) 



IV. EFFECTIVE SCATTERING AMPLITUDES 

We can apply the formulas derived in the previous sec- 
tion to the calculation of effective scattering times in both 
intrinsic and extrinsic semiconductors, provided they are 
in the degenerate regime (fcgT -C Ep). The intrinsic case 
will be considered first, since it is in this case that many- 
body effects play the dominant role. Electron-hole pairs 
are created in an intrinsic semiconductor by optical exci- 
tation and may, if the recombination time is sufficiently 
long, condense in a degenerate electron-hole liquid. The 
impurity concentration is negligible and electron-electron 
and, especially, electron-hole scattering plays the domi- 
nant role in controlling the DP mechanism of spin relax- 
ation. We will also consider, for completeness, the elec- 
tron liquid in extrinsic n-type doped semiconductors. In 
this case the doped electrons come from donor impurities 
and there is typically one impurity per electron. Further- 
more, because we are considering a bulk semiconductor, 
the impurities are homogeneously distributed throughout 
the electron liquid, and thus provide the dominant scat- 
tering mechanism - far more important than electron- 
electron scattering, as we will see. 

The crucial ingredient of the calculation is, of course, 
the effective interaction W to be used in Eqs. (54), (60), 
and (65). Here we have two options. The first option is to 
adopt a simple RPA-screened Coulomb interaction (also 
known as Lindhard screening) ; this leads to a parameter- 
free expression for the interaction, which should be exact 
in the high density limit. However the Lindhard screen- 
ing is known to be inaccurate as the density decreases, 
and furthermore the form of the RPA effective interac- 
tion misses important correlations between the electrons 
under consideration and the surrounding medium - cor- 
relations that become more and more important at low 
density. To counter these drawbacks one may resort to a 
second option, in which both the dielectric function and 
the effective interaction are modified by many-body local 
field factors which encapsulate exchange and correlation 
effects. (Discussion of these modified effective interac- 
tions can be found in Section 5.5 of reference [27].) Un- 
fortunately, even this approach is not problem-free, since 
the local field factors are imperfectly known in the elec- 
tron liquid, and even more so in the electron-hole liquid. 
Nevertheless, the behavior of the many-body local field 
factors is constrained by exact sum rules and limiting 
cases which, taken together, allow us to form a qualita- 
tively correct picture of the effective interaction at low 
density (still in the degenerate regime). We would like 
to point out that J. Zhou has also performed calcula- 



tions of the spin relaxation time including Singwi-Tosi- 
Land-Sjolander local field corrections, but in 2D GaAs 
systems. 28 

We define the two-particle bare Coulomb interaction 
(vij) and static Lindhard function (xoi) ahead of time: 



Vij(q) 



e r q' 

X 0i(q) = -Ni(0) 



1 1 

2 + ~ 4x 



x- 



ln 



1 - x 



(66) 
(67) 



where x = q/2k F , Ni(0) is the density of states at the 
Fermi level for electrons or holes, and e r is the relative 
dielectric constant of the medium. 

For the electron liquid, electron-electron interactions 
are handled by a spin-averaged Kukkoncn-Overhauser 
scattering amplitude: 29 

W^M = «n(«) + [vM (1 - GiMf Xii(q) , (68) 
with the static density-density response function: 

XoiO?) 



Xh(q) 



l-v n (q) [l-G!i(9)]Xoi(9) 



(69) 



Gfi(g) is the spin symmetric local field factor for elec- 
tron liquids. The scattering amplitude has the follow- 
ing physical interpretation: The first term vn(q) is the 
bare Coulomb interaction between electrons while the 
second term has an electron interacting with the den- 
sity induced in the electron liquid which in turn acts on 
another electron (each with a reduced Coulomb inter- 
action). Our calculations in the following section use a 
fit for G{ 1 (q) from Moroni, et. al. 30 (also reviewed in 
[31]) which is based on diffusive Monte Carlo studies of 
the static density-density response function. It should be 
noted that Gf 1 (g) used in these calculations is intended 
for use with 2 < r s < 10 and < q < 3k F , but lim- 
iting behavior and graphical analysis suggests the fit is 
still qualitatively reasonable for the range of r s exam- 
ined in this paper, 0.4 < r s < 2. Since we are restricted 
to the Fermi surface, the range of q is not a concern, 
< q < 2kp. Setting Gf 1 (g) = amounts to using the 
random phase approximation (RPA). 

Impurity scattering rates in the electron liquid are cal- 
culated with the electron-test charge effective scattering 
amplitude: 

W^(q) = v 21 (q) + «2i(g)xil («)«n (?) [1 - G' n (q)] , 

(70) 

where Xn(q) is as m Eq. (69) and we are assuming a 
uniform distribution of impurities. The physical inter- 
pretation is similar to that described for Eq. (68), except 
now the impurity interacts with the density induced in 
the electron liquid with the bare Coulomb interaction. 

For the electron-hole liquid, a multicomponent scat- 
tering amplitude analogous to the Kukkonen-Overhauser 
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formula is necessary: 
W% Hh (q)=v ij (q) 

2 

+ J2 v ik {q) [1 - G s ik (q)} Xkl{q)viM [l - Gfcfa)] 



kl=l 



(71) 



Elect ron- electron scattering is represented by VK™ L , and 
electron-hole scattering by W™ L . The density response 
function is found from 32 

[(x) _1 ]tf («) = [XOiiQT 1 Sij-ViM (1 - G'^q)) , (72) 

where x is the spin symmetric density response function 
in matrix form. The local field factors for electron-hole 
liquids (which should not be confused to be the same as 
those for the electron liquid) have not been well studied 
to the best of our knowledge. We opt to use an approxi- 
mate form of Gfj (q) from [32] with parameters based on 
a hole-to-electron mass ratio of mh/rn e — 6, which is ac- 
tually not too far from the value of the ratio in GaAs. 
Again, setting Gfj(q) = amounts to using RPA. 



V. CALCULATIONS OF THE EFFECTIVE 
SCATTERING RATES AND SPIN RELAXATION 
TIMES 

According to Eq. (26) to calculate the spin relaxation 
time we need two ingredients: the average scatter- 
ing time r* vg and the mean square of the spin-orbit field 
(&u)kr- The former has been calculated in Sec. Ill and 
is given by Eq. (54). The later can be found straightfor- 
wardly using the explicit form of f2k in Eq. (2): 



105 fr?E a 



(73) 



Inserting this result into Eq. (26) we get for the spin 
relaxation time 



1 

7W 



2^3 



32 qgg 

■m>h 2 E a 



(74) 



While the above results are valid for generic zincblcnde 
III-V semiconductors, we opt to include results for GaAs 
with the following properties: m c = 0.067m e , m v — 
0.47m e , E g = 1.43eV, e r = 13.2, and a c = 0.07. In 
Fig. 2, panel (a), the electron scattering rates from elec- 
trons and impurities in the electron liquid (EL) are plot- 
ted as functions of density. Panel (b) shows scatter- 
ing from electrons and holes in the electron-hole liquid 
(EHL). Recalling that the SRT is proportional to the 
scattering rate, we can see immediately via application 
of Matthiesen's rule that between e-e and e-i , the SRT 
will be controlled by e-i collisions in n-GaAs. Only at 
low densities do e-e collisions begin to make a significant 
contribution, and at that point degeneracy diminishes so 



that one should be more careful about including non- 
degenerate effects. In the electron-hole liquid (intrinsic 
GaAs), e-e collisions are downplayed even more and e- 
h collisions dominate the effective scattering rate. The 
physical reason for this is that holes, with their large mass 
and concurrently high density of states provide effective 
screening of the electron-electron interaction, which thus 
turns out to be much smaller than in the electron liq- 
uid. On the other hand, the presence of electron-hole 
scattering more than makes up for what is lost in the 
electron-electron scattering strength. 

When compared to the scattering rates of plane waves, 
the effective scattering rates are found to be generally 
smaller: see, again, Fig. 2. Approximating the effective 
scattering time by the plane wave lifetime can lead to a 
difference in the SRT of between and 1 order of mag- 
nitude. As discussed in the introduction, this happens 
because some collision processes cause a large change in 
momentum but a small change in the Dresselhaus field 
(consider for example a collision that takes us from a 
point in k-space where the Dresselhaus field vanishes to 
another symmetry-related point where it also vanishes); 
such a process contributes to the plane wave lifetime, but 
has no effect on spin relaxation. The use of the effective 
scattering rate rather than the plane wave approximation 
appears to be of utmost importance for high densities, 
but one must be cautious that the rate of spin relaxation 
does not approach the order of magnitude of the momen- 
tum scattering rate, as we will see momentarily. 



a) EL 



b) EHL 




15 20 25 30 



n (10 17 cirT 3 ) 



10 15 20 25 30 

n (10 17 cnT 3 ) 



FIG. 2. (Color online) The plane wave and effective scattering 
rates in the a) electron liquid (EL) and b) electron-hole liquid 
(EHL) have been calculated with RPA at T = 20K. 

In Fig. 3 the effect of local field factors on the effective 
scattering rates is apparent especially at low densities. 
While the local field factors for the electron-hole liquid 
are rough approximations based on m c /m v = 6, and so 
their results might not be quantitatively accurate, they 
likely reflect the qualitative trends as functions of den- 
sity. Additionally, while the local field factors for the elec- 
tron liquid have been calculated using modern quantum 
Monte Carlo analyses of the density response function, 
they were designed for the pure electron liquid (no im- 
purities) and considerably smaller densities. Disclaimers 
aside, we see that the inclusion of the local field fac- 
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tor generally enhances the effective scattering rates, and 
therefore suppresses the spin relaxation rate (e-e scatter- 
ing in the EHL is a bit of an anomaly here) . The enhance- 
ment of the interaction is expected on physical grounds 
since, as the density is lowered, the electron-hole liquid 
becomes increasingly "soft" - meaning a large density 
of low-lying excitations - and has a strong tendency to 
develop inhomogeneous density waves. The strong effec- 
tive interaction arising from this softening was recognized 
long ago as a possible mechanism of superconductivity in 
the electron-hole liquid. 32 However, the precise value of 
r s (the average inter-particle distance in units of the ef- 
fective Bohr radius) at which the transition would occur, 
as well as the actual degree of enhancement at a given 
r s , are still quite uncertain. 



a) EL b) EHL 




5 10 15 20 25 30 5 10 15 20 25 30 



n(10 17 cnT 3 ) n (10 17 cirT 3 ) 

FIG. 3. (Color online) a) Local field factors for the electron 
liquid (EL) are derived from fits by Moroni, et al. 30 and their 
inclusion in the scattering amplitude has the overall effect 
of enhancing the scattering rate, b) Local field factors for 
the electron-hole liquid (EHL) are based on a hole-to-electron 
mass ratio of 6 and are found in Vignale, et al. 32 . All rates 
have been calculated act T — 20K. 

With regard to Zhao's work mentioned in the intro- 
duction, the claims argued there are still relevant. The 
spin polarized packets in Zhao's experiment consist of a 
degenerate center and non-degenerate tails. In the tails, 
e-i collisions are reduced in liu of e-e collisions so that 
the effective scattering rate is generally smaller than in 
the center of the packet. The net effect is that spins relax 
faster in the tails than in the center. 

The spin relaxation times for n-GaAs and intrinsic 
GaAs are plotted in Fig. 4. As n-GaAs is largely dom- 
inated by e-i collisions, it can be understood why past 
theoretical curves fit the experimental data in references 
[33] and [34] so well. Caution should be taken when cal- 
culating the SRT for electrons in the electron-hole liq- 
uid. A region of validity becomes apparent for relaxation 
time approximation in Eq. (8). As pointed out earlier, 
SRTs shorter than the plane wave scattering time break 
the derivation of the effective scattering time. The plane 
wave lifetime in n-GaAs is sufficiently short to avoid this 
problem for all densities examined here, but e-h scat- 
tering has a comparatively long scattering lifetime which 
results in unfeasibly short SRTs. Admittedly, the over- 
all effective scattering rate ought to take into account 



all scatterers (phonons, impurities, etc.) which may ul- 
timately lengthen the SRT. The intrinsic GaAs example 
here is an idealized electron-hole liquid. 



a) n-GaAs b) Intrinsic GaAs 
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«(10 17 cm~ 3 ) n (10 17 cnT 3 ) 

FIG. 4. (Color online) The spin relaxation time in a) n- 
GaAs and b) intrinsic GaAs. Calculations were performed 
with RPA at T = 20K. n-GaAs is modeled by the EL and 
includes contributions from electron-electron and electron- 
impurity collisions. Impurity collisions dominate the SRT 
for the degenerate electron liquid. Intrinsic GaAs is modeled 
by the EHL and includes contributions from electron-electron 
and electron-hole collisions. Electron-hole collisions dominate 
the SRT for the degenerate electron-hole liquid. Where the 
SRT crosses the scattering time marks the beginning of the 
breakdown of the standard DP assumptions. From this point 
on, spin relaxation and momentum relaxation occur on simi- 
lar time scale and can no longer be separated. Eventually, at 
very high density, momentum relaxation becomes extremely 
slow and spin and momentum dynamics become decoupled 
again, this time with each spin performing an independent 
precession in the Dresselhaus field at a given point k. 



VI. CONCLUSION 

We have derived simple one- and two-dimensional in- 
tegrals for the effective scattering rate of electrons in the 
many body system, valid in the degenerate regime. The 
energy dependence has been handled according to ex- 
act expressions by Sykes/Brooker 20 and the angular de- 
pendence is a result of a Dresselhaus field acting as an 
addition spin-axis re-orientation mechanism. In general, 
scattering events with the Dresselhaus field present are 
less effective in randomizing the momentum axis, than 
without a field present. 

In highly degenerate systems, the contribution of 
electron-electron scattering to the spin relaxation time is 
minimal in comparison to electron-impurity or electron- 
hole contributions. For the case of high degeneracy in in- 
trinsic GaAs, we observe a limitation of the assumption 
that the timescale of spin relaxation is long compared to 
the scattering lifetime. In this regime the quasiparticles 
are essentially non-interacting and the DP mechanism is 
superseded by the spin precession of individual quasipar- 
ticles of essentially constant momentum. 

Local field factors, taking into account exchange 
and correlation effects, have been introduced into 
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the scattering rate calculations. Scattering rates are 
generally enhanced compared to using RPA, more-so at 
low densities, for both electron-electron interactions and 
electron-hole interactions. Significant improvements to 
the accuracy of the scattering rates could be made with 
new local field factors tailored to the densities of typical 
intrinsic and n-type III-V semiconductors. 



The addition theorem for spherical harmonics expands 
Pi in spherical harmonics: 

4-7T 

Pi (cos a) = — yr <p 9 ) yr ow*) ■ (A4) 

m 

Then, integration over dq is easy due to the orthogonality 
of spherical harmonics: 
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Appendix A: Equivalence of harmonics in Eq. (41) 
and the relaxation time approximation 



f 4tt 

=27rO k , z J d(cos a) WP 3 (cos a) . ( A5) 

That f2k,z can be extracted from each term in Eq. (41) 
goes to show that the relaxation time approximation in 
Eq. (7) and an angular-independent r£ are both reason- 
able claims. 



Appendix B: Reduction of the integral equation for 
the effective scattering time in electron-electron 
collisions 



It is relatively straightforward to demonstrate that in- 
tegration of 51 q 2 (where q can be any of k', p, or p') over 
e?q in Eq. (41) results in a term proportional to £1^^. Let 
us write ^ q ,z in terms of spherical harmonics: 

- fW^ll [Yi(# q ,<p q ) -Yf 2 (d q ,<p q )] (Al) 

where f2o is a constant. 

The angular portion of the collision integral in Eq. (41) 
has each f2 q ,z integrated with a scattering probability. 
These scattering probabilities can be expanded in Legen- 
dre polynomials of argument cos a = k • q: 

W =^2w t Pi (cos a) , (A2) 
I 

where W can represent either lW-kpk'p'1 2 OT 
|Wkpk'p'Wkp P 'k'|- The coefficients can later be 
found with 

21 + 1 f 1 

wi = — - — / d (cos a) WPi (cos a) . (A3) 
2 J-i 




FIG. 5. (Color online) For an arbitrary direction of k with 
respect to z, the direction of each momentum vector can be 
put entirely in terms of k and the angles 8, <fi, and tp p . Angle 
(j> defines the planes in which k,p and k',p' lie. 



We begin our solution for rj! by working from Eq. (41) 
with the replacements discussed, cubic symmetry allows 
Jlk x s — > £!k,z and angular independence in t£ allows 
T k "> T *k- 
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/ fe °(l - f° k )Q kiZ = E / dk ' d P - #)/ P (l fp') S ^ + £ p- ^ ~ V) ^k+p.k'+p' 

jl^kpk'p'l 2 [TfcOk,. - r fe *,O k ^] - i lWkpkvWkpp'k'1 t T fe^ z - rjf,n k /, z + r;ft p , z - J J . (Bl) 



The sum over p' can be evaluated immediately using 
conservation of momentum. Then, following a method 
introduced by Abrikosov and Khalatnikov 25 (or for an 
alternative derivation, Baym and Pethick 26 ), at low tem- 
peratures the integral over momentum space can be re- 
expressed as 



where 9 is the angle between k and p, and is a polar 
angle about the k + p axis (See Fig. 5). Aside from spin 
dependence, the scattering amplitudes are functions only 
of 9 and 6. 



dk' dp 



2h 6 cos (9/2) ^ £k ' d£p ^ £p ' Sin ^ ^ ^ p ' 

" (B2) 



At this point our kinetic equation has the form 







2ml 



k,z 



W(2ixY 



de k > 



ds r 



de r . 



" % ) - fp>)S (£fc + £ P - £fc' - V) / d<jy I d9 I d<p p 



sin(0/2) j|Wk pkV | 2 [O kjZ T fc * - n k ',zTfe'] - ^ |^ kpk 'p'T^ k pp' k '| [^k^ - f^,^*' + ^ p , z T p * - fi p ', z T*,] | . 



(B3) 



The spherical harmonic definition of f2 kz from Eq. (42) 
is useful at this point. The orthogonality associated with 
the spherical harmonic functions is exploited by multi- 
plying through by £l kz in Eq. (B3) and integrating over 
siwdkdflk difik- The left hand side reduces with use of 

J dd k J d^ k smi9 k Ql z = C 2 k e —, (B4) 

where C 2 = a^ 4 /(2m^ 9 ). 

There is still the matter of putting the angles corre- 
sponding to vectors p, k', and p' in in terms of the inte- 
gration variables. This is accomplished by introducing a 
right hand orthonormal basis of vectors: 

k = {sin# fe cos<pfc,sin# fe sin<pk,cosi?k} , (B5a) 
t = {cosz9fc cos<y9fc, cos'&k sin ip k, — sin^fc} , (B5b) 
n = k x t . (B5c) 

Aligning k along z, it is seen that 

p = cos 9k + sin 9 cos (p p t + sin 9 sin ip p ii . (B6) 

A rotation formula can be applied so that vectors k and 



I 

p are turned about the k + p axis to find k' and p': 
t^ + ^Wg^)**, (B7> 

Standard relations are used to put the angular compo- 
nents of k', p, and p' in terms of the integration variables: 

cos -d q = q z , (B9) 

cos ^ = ^=^- ( B1 °) 

This concludes the needed angular transformations. 
Let us define the angular-only portion of fi kz with 
K("& k ,(p k ). In other words 

K(i?k,<£k) = cos d k sin 2 fi k cos(2ip k ) 

= (cosi? fe -cos 3 i? fc )(2cos 2 ^ fe - 1) , (Bll) 

where identities have been used to write k in terms of 
cosine functions for direct use of Eqs. (B9) and (BIO). 
This leaves us with 
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16tt 
105 



2mt 



fi 7 (27r) 5 



de k > 



d£r. 



' " f %)f P {l-f p .)He k +e p -E k ,~e p ,) 



Jk 



j dip k / d&k / d(j> d6 d(p p smi9 k sm(9/2)K(-d k ,tp k ) i | W k k' 
Jo Jo Jo Jo Jo I 



K{d k ,ip k )rl 



{k/k 1 



,3 T fe' 



- g |WWW kp , 



K(6 k ,ip k )T k - 3 Ty 



(fc/p) 3 ^ 



(fc/p') 3 



(B12) 



Here, we have made use of the fact that the scattering 
amplitude in the direct portion of the collision integral 
only depends on the momentum transfer k — > k', whereas 
the exchange portion includes the possibilities of k — > k' 
and k — > p'. For scattering near the Fermi surface, the 
scattering amplitudes only vary with the angle between 
incoming and outgoing momenta. These angles can be 
written in terms of just 9 and <j) with the following rela- 



tions: 

cos #i = k • k' 

1 (B13) 
= -(1 + cos 9 + cos <j) — cos 9 cos 0), 

cos #2 = k • p' 

1 (B14) 
= -(1 + cos 9 — cos + cos 9 cos 0). 

The scattering amplitudes written as functions of 9\ and 
92 are 

WW=W(cos0i), (B15) 
=W(cos0 2 ). (B16) 

For reference, the dimensionless momentum transfers are 

gi = |k - k'|/ (2k F ) = sin (0/2) sin (0/2) , (B17) 
<Z2 = |k - p'|/ (2fc F ) = sin (6>/2) |cos (<j>/2)\ . (B18) 

The absolute value is a result of the range of integration, 
< <t> < 2tt. 

Integration over dip p , d-d k and d(p k can be done imme- 
diately with the result 



h 7 (2ir) 4 



mi 



/OO POO POO / 1 — f \ 

de k , \ de p \ de p , 
-oo J — oo J — oo V 1 ^ fc / 



IW, 



kk' 



cos(0/2) 



P 3 (cO8 0i) 



(jfe/jf) 



A3 



2 cos(0/2) 



P 3 (cosgi) , 

(fc/fc') 3 TA 



- ey - e p -) y 



P 3 (cos0) t P 3 (cos0 2 ) 
(fc/p) 3 Tp ~ (fc/j/) 3 



, (B19) 



where = sin d6d<p should not be confused with the 
Dresselhaus Larmor frequency. 

Notice that without the Legendre polynomial terms, 
the scattering time could be factored out of the inte- 
gral and we would have exactly the scattering rate of 
plane waves, as can be referenced in [27]. Given that 
well known solution, we are led to believe t£ will have 
even dependences on energy and temperature. 

The ratios (k/k 1 ) 3 , (k/p) 3 , and (k/p') 3 vary slowly 
across the Fermi surface when compared to the Fermi- 



Dirac functions and anticipated e k — \i dependence in t£ . 
They are all set to 1. We shall consider only static screen- 
ing in the scattering amplitudes, a reasonable assumption 
at low temperature. The resulting integral equation for 
has been exactly solved by Sykes and Brooker. 20 We 
will reproduce a simplified version of the solution here. 



The energy integrals are written in terms of unitless 
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variables: 



and the <5-function evaluated for (e p — /j,)/k B T to give 



£ = (e fc -/i)/fc B T, 
x = (e k , - fi)/k B T, 
y = (v - v)/k B T , 



(B20) 
(B21) 
(B22) 



n 7 (27r) 4 
ml (k B Tf 



dx 



dy 



/"(-*) 

/°(-0 



f°(x + y-0.f°(-y) / dQ 



kk' 



cos(6»/2) 



[T*(O-P3(cos0i)T*(aO] 



"2 cos(0/2) 



[r* (0 - P 3 (cos ^) r* (a;) + P 3 (cos 0) r* (z + y - - P 3 (cos 9 2 ) T 



*(»)]} • 



(B23) 



The Fermi-Dirac functions are now described by 

A*) ' 



1 



(B24) 



Integration of t*(x + y — £) over dy is equivalent to in- 
tegration of t*(— y) over dy. This relation, along with 
a simple swap of x <-> y variables in terms which have 
r*(±y), allows the integration over dy to be performed: 



J — ( 



dy/ (x + y-e)/°(-y) = 



1 - e 



(B25) 



With an even dependence on x, we can also let r* (—or) — > 
r*(x). To simplify the form of the integral equation, the 
angular integrals are represented by constants B and A, 
and the remaining x-dependent factor by K(x,£): 



A 1 
A 2 
Ai 
A 2 

A 



f°i-x} r g-x 

7°(-0 |.l-e*-« 



dfi 

dfi 



IW, 



kk' 



/I 
1 



cos(0/2) ' 

cos(0/2) ' 

dfi-^^P 3 (cos0i) , 
cos(#/2) 

cos(fl/2) [P3(C ° Sgl) 



-P 3 (cos0)+P 3 (cos0 2 )] , 

m 3 c (k B T) V 2 
A1A1- A 2 A 2 /2 



Ai - A 2 /2 



(B26) 
(B27) 
(B28) 
(B29) 

(B30) 

(B31) 
(B32) 



The remaining integral equation for t*(£) is exactly of 
the form in Sykes and Brooker: 

/oo 
dx [t* (0 - At* (a;)] . (B33) 
-00 



Notice that the kernel of this integral equation involves 
only two-dimensional angular integrals on the Fermi sur- 
face. In fact, since A\ and Ai do not contain any ex- 
change terms, they can actually be simplified to one- 
dimensional integrals of the momentum transfer q\ in 
Eq. (B17): 



! =8tt / dqi 
Jo 

/ dqx 

Jo 



A 

Ai =8tt 



|W(?i)| 2 , (B34) 



\W(q i rP 3 (l-2q()). (B35) 



The solution for r*(£) is found by converting the inte- 
gral equation to a differential equation via Fourier trans- 
form, utilizing the convolution theorem. For r*(£) with 
an even dependence on £, 



where 



cosh(£/2) 
2^ 



-00 in 



(4i + 3)$2tM 

-A) ' 

(B36) 



" ^ A 2 ; (A 2( 



=P|+i (tanhTrw) , 
A, = £(i + l)(i + 2), 



(B37) 
(B38) 



and p™(x) are the associated Legendre polynomials. 
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